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■ ABSTRACT 

' We present a detailed study of the spatially resolved kinematics, star-formation and stellar 

, mass in a highly amplified galaxy at z = 4.92 behind the lensing cluster MS 1358+62. We 

use the observed optical, near- and mid-infrared imaging from HST ACS & NICMOS and 
Spitzer IRAC to derive the stellar mass and the Gemini/NIFS IFU to investigate the veloc- 
ity structure of the galaxy from the nebular [Oll]AA3726.8,3728.9 emission. Using a detailed 
gravitational lens model, we account for lensing amplification factor 12.5±2.0 and find that 
this intrinsically L* galaxy has a stellar mass of M*=7±2x 1O^M0, a dynamical mass of 
M<jy„=3±lxlO^c.sc2('/jM0 (within of 2kpc) and a star-formation rate of 42±8M0yr"^ 
The source-plane UV/optical morphology of this galaxy is dominated by five discrete star- 
forming regions. Exploiting the dynamical information we derive masses for individual star- 

■ forming regions of Md '^IO^^^Mq with sizes of ~200pc. We find that, at a fixed size, the 
I Star-formation rate density within these Hll regions is approximately two orders of magni- 
. tude greater than those observed in local spiral/starburst galaxies, but consistent with the most 

massive Hll regions in the local Universe such as 30Doradus. Finally, we compare the spa- 
tially resolved nebular emission line velocity with the Lya and UV-ISM lines and find that 
. this galaxy is surrounded by a galactic scale outflow in which the Lya appears redshifted 

^ ' by ^150kms^^ and the UV-ISM lines blue-shifted by ^-200kms^^ from the (systemic) 

\ nebular emission. The velocity structure of the outflow mirrors that of the nebular emission 

suggesting the outflow is young (< 15 Myr), and has yet to burst out of the system. Taken to- 
gether, these results suggest that this young galaxy is undergoing its first major epoch of mass 
assembly. 

Key words: galaxies: high-redshift — galaxies: evolution; — galaxies: star formation rates, 
dynamics, — galaxies: individual: MS 1358+62arc 



1 INTRODUCTION 

Deep imaging surveys, particularly with Hubble Space Telescope 
(HST) have now uncovered tho usands of star-f orming g alaxies in 
the redshift range z = 4 — 6 (eg. Ouchi et al.l200 4: Giavali sco et al 
20041 : lYoshida et ai] |2006| ; iBeckwith et al.[|2006l : iBouwens et al 
2007h . Studying such galaxies are key steps towards understand- 
ing the physical processes that drive galaxy formation in the early 
Universe, probing the properties of galaxies when they formed 
their first generation of stars. Concentrated follow-up of Lyman- 
break galaxies at this epoch has provided estimates of stellar 



masses (^lO^-i^Mo), ages (~150-300Myr) and hence inferred 
star-formation ra t es of a few lO's yr~ ^ jVerma et al.l 1 20071 : 
IStark e tal.' 20091; iMcLure et all '2009). With high specific star- 
formation rates, these young galaxies are likely assembling their 
first significant stellar mass. Moreover, since this epoch is within 
~300Myr of the end of reionisation (assuming reionisation was 
completed by z ~6) many of these galaxies should have been ac- 
tively star-forming at z ^7. As such, probing the ubiquity of star- 
forming galaxies, the ionising photon density from star-formation 
and their relation to AGN are important issues, potentially allow- 
ing us to understand how and when the first galaxies were assem- 
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bled and how the Universe was reionised. Moreover, theoretical 
models suggest that the most rapid epoch of dark h alo assembly 
for g alaxies such as the Milky- Way occur at z ^5 ( IMo & White! 
12002 : Okamoto et al. 2005). Through detailed studies of the phys- 
ical properties of galaxies at this epoch (such as masses and star- 
formation rates), we can accurately assess their contribution to the 
mass growth of galaxies like the Milky- Way, and their effect on the 
gaseous IGM. 

Indeed, one of the key findings in recent years is the dis- 
covery that a significant fraction of high-red shift galaxies are 
surrounded by "superwind s" (e.g. Franx et al. '1997; Pettini et al.' 
2002b; Shapley et al. 20031; lBower"e t al. 2004; Wilman et al. 2005; 
Swinbank et al.l 120071) - starburst and/or AGN driven outflows 



which expel gas from the galaxy potential and hence play no fur- 
ther role in the star-formation history of the galaxy. Such feed- 
back processes may offer natural explanations to the shape and 
normalisation of the local luminosity function (Benson et al. 2003; 
iBower et alj|2006l ; ICroton et al.ll2006l) and offer natural explana- 
tions as to why only 10% of the bary ons cool to form stars 
jWhite & Ree3ll978l ; Is^ogh et al.ll200ll) . Although velocity off- 
sets between star-forming regions and outflowing gas of several 
hundred km/s have been measured (suggestive of large scale out- 
flows), constraining the geometry and mass loading is vital to test 
whether the outflows are truly large-scale (as often obse rved in 
low-r edshift Ultra-Luminous Infrared Galaxies (ULIRGS); lMartiril 
|200^) . or whether the outflows are confined to individual star- 
forming regions. Such constraints are vital in order to test whether 
the outflowing material escapes into the IGM, or whether it even- 
tually stalls, fragments and drains back down onto the galaxy disk, 
causing further bursts of star-formation. 

However, probing the masses, star-formation rates, and in- 
teraction between sta r-formation a nd gas dyna mics which are 
routine at z ^2 (eg. Genzel et al. 1 2006; Fors ter Schreiber et al.l 
2006'; 'Swinbank et al.' 2006b; Law etal J 120071; Istark et alj|200i ; 
Lehnert et al. 2009; Law et al. 2009) are difficult beyond z ~3 due 
to a combination of surface brightness dimming and the smaller 
physical sizes of the earliest systems (and the strong evolution 
in the UV luminosity functi on between z ^ 3-6 means there are 
fewer luminous sources; eg. 's tark et alj |200^. Prior to the launch 
of James Webb Space Telescope (JWST) one particularly appeal- 
ing route to study the internal structures of primitive galaxies at 
z ~ 5 is to target galaxies which have been highly magnified 
by their serendipitous alignment with a foreground massive clus- 
ter. The lensing magnification (up to a factor 30 x in the most 
extreme cases) has two effects: the galaxy image is amplified at 
a fixed surface brightness so that the galaxy image appears both 
brighter and larger. By combining gravitational lensing with re- 
solved spectroscopic imaging, detailed maps of the sizes, dynamics 
and hence masses of both young galaxies and even starburst/HlI re- 
gions within these galaxies can be made. Such observations can 
be used to investigate whether the physical conditions for star- 
formation were significantly different from comparably luminous 
galaxies at z ~ 2 and galaxies today. Indeed, gravitational lens- 
ing and resolved spectroscopic imaging has been used to inves- 
tigate galaxy dynamics on ^Ikpc scales from z = 1-3 (eg. 
Swinbank et ai]|2003l. l2006al ; ISalucci et al.]|2007l : iNesvadba etatl 
2006 1). whilst coupling lensing with adaptive optics (AO) source- 
plane studies can even reach <100pc (Stark et al. 2008; Jones et 
al. 2009 in prep). 

Studying the physical processes of galaxy formation on scales 
comparable to the largest Hll regions is likely to have broad impli- 
cations for our understanding of the dominant mechanism by which 



galaxies at high redshift assemble their stellar mass, and what 
proce sses drive the star-formation activity dBournaud & Elmegreenl 
I2OO9'). For example, recent hydro-dynamical simulations have sug- 
gested that continuous gas accretion (so-called "cold-flows" ) may 
play an important role in driving star-formation at 2 > 2 jPavd 
l2008l ; lDekeletalJl200 9'). In such models, these cold streams keep 
the rotating disk configuration intact, but with turbulent star- 
forming clumps which eventually merge to form a central spheroid. 
This is a provocative result, as it has generally been accepted that 
major mergers are the dominant m echanism by which m assive 
galaxies assemble their stellar mass dConselice et alj|2008l) . How- 
ever, to test such models in detail, the galaxy velocity field and the 
turbulent speeds on scales comparable to individual Hll regions in 
a large sample of high-redshift, star-forming galaxies are required. 
Resolved dynamics of galaxies at 2 ~ 2 have lent support to the 
cold-stream model, but at the spatial scales available in non-lensed 
galaxies (I~4 kpc) it is difficult to draw definitive conclusions. Nev- 
ertheless, the advances in instrumentation (particularly sensitive 
integral-field units and AO) are beginning to allow studies of the 
kinematics, distribution of star-formation, gas and stars and chem- 
ical properties in high-redshift star-forming galaxies which can be 
fed into (and thus test) numerical simulations. 

In this paper we report the detailed follow-up of a lensed 
galaxy at z = 4.92 be hind the len s ing c luster MS 1358-1-62, 
first studied in detail by iFranx et al] ( Il997l) . We use the HST 
ACS & NICMOS imaging to investigate the rest-frame UV/optical 
morphology. Together with Spitzer IRAC imaging we construct 
the spectral energy distribution (SED) to estimate the stellar 
mass. Using the Gemini/NIFS integral field unit we map the 
[Oll]AA3726.2,3728.9 emission line strength and dynamics, distri- 
bution of star-formation and outflow energetics. We have adopted 
a standai'd ACDM cosmology with Ho = 72 km s"\ = 0.27 
and Qa ~ 0.73. All quoted magnitudes are on the AB system. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Optical", near- and mid-infrared imaging 

Details of the optical, ne ar- and mid- infrared observations of this 
galaxy cluster are given in lRichard et alj j2008l) . Briefly, the lensing 
cluster MS 1358-1-62 (Zrf=0.33; a: 13:59:50 5: -1-62:31:05, J2000) 
was observed with HST ACS and NICMOS as part of program 
PID: 9717 & 10504 using the 435W, 475W, 625W, 775W, 850LP 
FllOW and F160W filters with integration times 5.4-1 6 ks each 
in the optical bands and ~I9ks each in the near-infrared bands. 
The optical and near-infrared imaging reach typical 5a depths 
of AB=26 .5-27.5 magnitudes . We also include if-band photom- 
etry from ISoifer et"al1 ( Il998[) which reaches a 5(j sensitivity of 
Kab ~23. The Spitzer IRAC 3.6, 4.5, 5.8 and 8.0^m imaging 
comprises 2.4 ks in each band, reaching 5a depths of AB~24 
dRichard et"ai]|2008l) . 

The z = 4.92 arc lies approximately 21" south west of the 
brightest cluster galaxy at a: 13:59:48.7 5: -1-62:30:48.34 (J2000) 
and was discovere d during a spectroscopic survey of this cluster by 
iFranx et aO d 1997b in which the redshifted Lya and UV-ISM lines 
were used to derive an unambiguous redshift of 2 = 4.92 (a com- 
panion galaxy at the same redshift located approximately 200 kpc 
away in the source-plane was also discovered during the same sur- 
vey). In Fig. [T] we show a true-colour HST ACS V I z—hand im- 
age of the cluster core around the 2 = 4.92 arc. We also show 
an ACS/NICMOS / J_ff-band true-colour image, as well as thumb- 
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Figure 1. Trae colour HST ACS y_R/-band image of the lensed z = 4.92 galaxy behind MS 1358+62. The image is centered at the position of a bright cluster 
elliptical at a =13:39:49.51 <5 =+62:30:48.69 (note that the BCG lies approximately 20" south). The z = 4.92 critical curve from our best-fit lens model 
is also overlaid. The HST imaging clearly shows the galaxy resolved into two images of the background galaxy, with the morphology dominated by up to six 
star-forming regions surrounded by a diffuse halo. The dashed white box shows the NIFS IFU field of view used to cover the z = 4.92 galaxy image. The 
upper left insets show a colour image generated from the NICMOS FllO and F160 images showing that the J- and /f-band morphologies clearly trace the 
observed yi?/-band morphology. The upper middle panel shows the white light (wavelength collapsed) image of the [Oil] emission from the IFU datacube 
(dark grey denotes highest intensity). The boxes (labeled 0-4) denote the regions from which spectra were extracted (see Fig. 4). The upper right insets show 
the HST ACS /-band image (dark grey denotes high intensity) with contours from the [Oil] emission line image overlaid, showing that the [Oil] emission line 
moiphology follows that seen in the rest-frame UV. The contours denote [Oil] emission line flux levels of 0.2-1 X lO"^'^ erg s"'^ cm~^/sq-arcsec in increasing 
in units of 0.1 X 10^ erg s^^ cm~^/sq-arcsec. The lower four insets show the Spitzer IRAC 3.6, 4.5, 5.8 and 8.0/^m images around the lensed galaxy. In each 
panel we mark the position of the brightest component of the ^ = 4.92 arc, which is clearly visible in all four channels (each of these panels is 30" on a side). 



4 



nails around the galaxy in the four IRAC images showing that the 
galaxy is detected from l/-band to the IRAC 8.0/im channel. 

To calculate the photometry of the lensed galaxy we first con- 
sider the bright elliptical galaxy locate ~6" due East (which has 
a low surface brightness halo which extends beyond the z — 4.92 
galaxy image). We construct a model of the elliptical galaxy sur- 
face photometry using th e IRAF BMODEL ellipse fitting algorithm 
jjedrzeiewski et aUll987l) (we note that we mask the z = 4.92 
galaxy image during the fit). In this model, the surface photome- 
try includes the C4 Fourier coefficient (which is necessary to de- 
scribe the boxiness of the isophotes), but all other Fourier terms o 
were forced to zero. We then subtract the best fit model from the ^ 
galaxy image and then use SEXTRACTOR dSertin & Amouts|[l993) 
to estimate the residual background within the frame. Since accu- 
rate colours are required to derive the spectral energy distribution 
we calculate the magnitude of the lensed galaxy in the various pass 
bands using an elliptical aperture (approximately 3"xl.5" cen- 
tered on the first image of the arc as covered by the NIFS IFU). 
This aperture is constructed from the coadded /- and z-band im- 
ages since these represent the highest signal-to-noise in the z=4.92 
arc. The same aperture is then applied to all of the images to mea- 
sure the photometry. In the IRAC channels we apply an aperture 
correction to account for the flux lying outside the aperture (~0.6- 
0.8 magnitudes from 3.6-8.0/im) as determined from bright (unsat- 
urated) point sources in the image. 

In Table I we give the observed photometry of the galaxy from 
the HST ACS & NICMOS and Spitzer IRAC imaging and in Fig. [3] 
we show the observed photometry with the best fit model SED over- 
laid (see ii l3.2t . We note that the large error-bars (particularly at 5.8 
and 8.0/im) reflect the uncertainty in recovering the photometry 
from the arc after modeling and removing the foreground elliptical 
galaxy which dominates the light, particularly at the reddest wave- 
lengths. We also note that there is MIPS 24jim coverage of this 
cluster with a bright source at the position of the galaxy. However, 
due to the blending with the foreground lensing galaxy (located 
~6" East), it is not possible to reliably deblend the photometry 
from the elliptical galaxy. 

2.2 Spectroscopic Imaging 

Three dimensional spectroscopic imaging observations around 
the redshifted [Oll]AA3726.2,3728.9A emission line were taken 
with the Gemini-North Near-Infrared Integral Field Spectrometer 
(NIFS) between 2006 February 08 and 2006 February 09 during 
science verification time. 

The Gemini-NIFS IFU uses an image sheer to take a 
3.0" X 3.0" field and divides it into 29 slices of width 0.103". The 
dispersed spectra from the slices are reformatted on the detector 
to provide 2-dimensional spectro-imaging, in our case using the 
K-band grism covering a wavelength range of 2.00-2.43 /xm. The 
observations were performed using an ABC configuration in which 
we chopped by 6" to blank sky to achieve sky subtraction. Indi- 
vidual exposures were 600 s and each observing block was 3.6 ks, 
which was repeated four times for resulting in an integration time 
of 14.4 ks (with a total of 9.6 ks on source and 4.8 ks on sky). 

We reduced the data with the standard Gemini IRAF NIFS 
pipeline which includes extraction, sky subtraction, wavelength 
calibration, and flat-fielding. Residual OH sky emission lin es were 
remov ed using the sky-subtraction techniques described in lOaviej 
( I2OO7I) . To accurately align and mosaic the individual datacubes 
we created white light (wavelength collapsed) images around the 
redshifted [Oil] line from each observing block and centroid the 
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Figure 2. True colour HST ACS VRI-band reconstruction of the lensed 
galaxy using the mass model described in §[5] The amplification of the 
galaxy is a factor 12.5±2.0x. In the source plane, the galaxy has a spatial 
extent of ~4kpc and comprises at least five discrete star-forming clumps. 
The largest of these is only marginally resolved, with a rest-frame optical 
half fight radii of ~200 pc. 

galaxy within the data cube. These were then spatially aligned and 
co-added using an average with a 3a clipping threshold to remove 
remaining cosmetic defects and cosmic rays. Flux calibration was 
carried out by observing a bright AOV standard star (HIP 78017) at 
similar airmass to the target galaxy immediately after each observ- 
ing block. From the reduced standard star cube, we measure a K- 
band seeing of FWHM=0.45". We note that accounting for lensing 
amplification factor I2.5±2.0x (derived in ij3} this corresponds to 
a source-plane FWHM of ~0.6kpc. The spectral resolution of the 
data (measured from the sky-lines at ~2.2/im) is _R = A/AA=5300 
which corresponds to (7=3 A or 25 kms^^. In all following sections, 
quoted line widths are deconvolved for the instrumental resolution. 



3 ANALYSIS 
3.1 Lens modeling 

We first estimate the amplification of the galaxy by making use 
of the updated gravitational lens model from lRichard et al. (200^. 
Fig. [T] shows the cluster region around the lensed galaxy with 
the best-fit z = 4.92 critical curve overlaid. Using the map- 
ping between image and source-plane co-ordinates, we ray-trace 
the 2 = 4.92 galaxy image and reconstruct the source-plane 
morphology, deriving a luminosity weighted magnification factor 
/i=I2.5±2.0x (which corresponds to Am=2.7±0.2 magnitudes). 
The error-bar on /i is derived by ray-tracing the family of ac- 
ceptable gravitational lens models which adequately reproduce the 
multiply-imaged galaxies in the cluster. 
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Figure 3. Observed Spectral Energy Distribution (SED) of the galaxy from 
the observed HST ACS/NICMOS and Spiner IRAC photometry. We over- 
lay the best-fit SED from H YPERZ which predicts a young age of ~20 Myr 
and moderate (Ai,=0.25) reddening. Fitting the SED using we derive a stel- 
lar mass (corrected for lensing amplification) of = 7 ± 2 X 10* Mq. 



Accounting for lensing amplification, the intrinsic magnitude 
of the galaxy is Iab = 24.94 ± 0.08. For comparison, deep imag- 
ing surveys of Lyman-break galaxies at 2: ^-^5 have derived an L* 
of lAB ~ 25.3 suggesting that the z — 4.92 galaxy studied here is 
typical of th e UV-continuum selected samples at this redshift (eg. 
lOuch i et al. 20041. 

In Fig.[2]we show the reconstructed source-plane image of the 
galaxy which shows that the rest-frame UV/optical morphology is 
elongated, but is dominated by two star-forming regions separated 
by ~4 kpc with at least four more lower-surface brightness regions 
within the galaxy and surrounded by a low surface brightness halo. 

To estimate the size of the source-plane PSF in each band, 
we first measure the size of the PSF using the i? ~19 star located 
5" due south of the arc, which has a FWHM=0.10, 0.33, 0.30" 
in the observed I-, J- and iJ-bands respectively. The two largest 
star-forming regions (at approximately [1.8,2.0] and [-1.0,-2.1] kpc 
in Fig. O have an /-band FWHM which is marginally resolved 
(FWHM=0.17±0.03"), but both are unresolved in J and -bands. 
Since the lensing amplification of the source has a preferential di- 
rection, we reconstruct the source-plane image of the star at the lo- 
cations of the brightest star-forming regions and find that the ACS 
PSF transforms to 0.06" x0.014", (320x85 pc). Both of the bright- 
est two Hll regions are marginally resolved in the /-band imaging 
(FWHM=0.17±0.03" and FWHM=0.15±0.04 for regions A and 
B in Fig. [T] respectively) corresponding to FWHM=500±100pc 
and FWHM=400±150pc (deconvolved for PSF). We note that the 
size of the largest of these two Hll regions is comparable to the 
size measured by Franx et al. ( 1997) (FWHM~390pc) from HST 
WFPC2 observations. In the J- and //-bands, the source-plane 
PSFs is approximately 800x220pc and 300x1 kpc in the NIFS 
[On] emission line map. We therefore conservatively suggests sizes 
of FWHM<500 and 400 pc in J- and //-bands respectively and 
<600 pc in the On emission line map for the two Hit regions, A & 
B respectively. We discuss the implications of these sizes in §3.6, 
but note that clearly higher spatial resolution IFU observations us- 
ing laser-guide star AO would be very valuable to further constrain 
their sizes. 
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Figure 4. One-dimensional spectrum of the five star-forming regions within 
the 2 = 4.92 galaxy from the NIFS IFU observations. In all panels we 
show the position of the [OlI]AA3726.8,3728.9 doublet at a fixed redshift of 
z =4.9296. In each panel we also overlay the best fitting spectrum and give 
the (rest-frame) velocity offset from the systemic redshift and line width of 
the best fitting double Gaussian profile. From the spectra, the velocity gra- 
dient of Av=180±20 km can easily be seen. The final panel shows the 
one-dime nsional spectrum o f the 2 = 4.92 galaxy around the Lya emission 
line ffrom lFranx et al.l 19971) showing the asymmetric line profile. Using the 
systemic redshift measured from the nebular emission, we de-redshifted the 
spectrum and find that the asymmetric emission line is best fit with an un- 
derlying Gaussian emission line profile (blue dashed line) combined with a 
Voigt absorber in the blue wing. 



3.2 Stellar mass 



ISoifer et"al] ( 1 19981) use optical and near-infrared photometry to es- 
timate the stellar mass in this galaxy, deriving M*=2-16x IO^Mq. 
Here, we make use of the much deeper ACS/NICMOS and, partic- 
ularly IRAC imaging in order to improve the estimate of the stellar 
mass. We infer the stellar mass of the galaxy by fitting the latest 
Chariot & Bruzual (2007) (hereafter CB07) stellar population syn- 
thesis models to the observed SED (which include an improved 
treatment of the thermally-pulsating asymptotic giant branch (TP- 
AGB) phase of stellar evolution). Following Stark et al. ( 2009), we 
allow 200 age steps f rom 10^ t o 10^ y r (approximately logarithmi- 
cally spaced), with a lSalpeteil ( fT953) initial mass function (IMF) 
and consider both a solar and sub-solar (0.2 Zq) metallicity. We fit 
both constant star-formation histories and exponentially decaying 
(r) models with e-folding decay times of 100 Myr. We consider the 
effects of dust in the modeling by adopting the lCalzetti et alj ( [l994h 
reddening law (using colour excesses of E{B — V)=0.0, 0.03, 0.05, 
0.1, 0.3 & 1.0). The resulting template SEDs are convolved with the 
filter response curves, and by comparing the observed magnitudes 
and their associated errors, we compute a reduced value for each 
set of parameters (age, E{B — V) and normalistation). The la un- 
certainties are calculated from the parameter range which produce 
Ax^=l with respect to the best fit SED (whilst marginalising over 
the other free parameters). 

For consi stency, we fi r st mod el the IJH K-band galaxy pho- 
tometry from ISoifer et"aT I ( Il998l) . deriving acceptable solutions 



with E{B 
of M 



y)=0.06-0.4, ages= 1 X 10®"^ Myr and stellar masses 



1 X lO^^x 10^" Mq (which is consistent with Soi fer et all 
1 19981) . However, the much deeper //5r ACS/NICMOS and Spitzer 
photometry (VIzJHK & IRAC 3. 6,4.5,5. 8,8.0Aim) allows much 
better constraints to be derived, and we infer a best fit SED with an 
age of 14±7 Myr, an E{B - V)=0.05±0.05 and a stellar mass of 



6 



M*=7±2x 10** M0, with a 0.2 metallicity template marginally 
prefered over solar metallicity (although formally the difference in 
Ax'^ is less than 1 a between the two). We note that the constant 
star-formation and r models provide equally good fits to the data, 
with stellar masses which vary by <5%. We also attempt fits ex- 
cluding the IRAC 5.8 & S.O/^m photometry, but the resulting stellar 
mass and age agree within the \o errors. We note that the stellar 
mass derived here is significantly lower (at least a factor >3 x) than 
that derived by Soifer et al. 1998 , which most likely arises due to 
the improved photometric constraints and very young age as sug- 
gested by the IRAC photometry. 

3.3 Star-formation rate and dynamics 

Next we use the spatially resolved [Oll]AA3726.8,3728.9 emis- 
sion line maps from our NIFS IFU observation to estimate the 
distribution of star-formation and dynamics within the galaxy. 
As Fig. [T] shows, the nebular line emission is clearly patchy, 
and dominated by the brightest regions seen in the HST 
ACS and NICMOS imaging. Collapsing the brightest regions 
within the galaxy we find that the integrated [Oil] emis- 
sion line flux is /[o//] = 1.5±0.3 x 10"^'' ergcm~^s^^. Ac- 
counting for an amplification factor (^=12.5±2.0x) and us- 
ing the conversion b etween [O il ] emi ssion line flux and star- 
formation rate from iKennicutll ( Il998h we derive an intrinsic 
star-formation rate of 5'Fi?[o//]=42±8 Mq yr~^ (uncorrected for 
reddening). This star-formation rate is consistent with the 3(t 
upper-limit of SFR~80 M(:) yr~^ (corrected f or lensing) from 
SCUBA 850^m imaging jKnudsen et aLll2008 h. We note that if 
the star-formation rate derived from the [Oil] emission line flux 
(S'F_R=42±8 Mq yr~^) has been sustained, then it takes just 
~15Myr to build a stellar mass of M*=7x 10** Mq, suggesting 
the current burst may be the first major episode of star-formation 
within this galaxy (which is also consistent with the young SED). 
It is also useful to note that the starburst intensity (SFR/area) for 
the galaxy (Se ~4x 10^^ L© kpc~'^) is com parable with local- (eg . 
Meurer et al]ll997h and high redshift- (eg. loiavaUsco et al.|[l996l : 
Steidel et aljl996l) galaxies, suggesting the same mechanisms may 
limit the starburst intensity at z ~5 as those at lower redshifts. 

We also use the spatially resolved [Oil] emission line doublet 
from the NIFS IFU observations to investigate the dynamics of the 
galaxy. By collapsing the velocity field from the brightest regions 
within the datacube we derive a peak-to-peak velocity gradient of 
Aw=180±20km s~^ across 4 kpc in projection (Fig.O. If these dy- 
namics reflect virialised motion within a bound system, then we es- 
timate a dynamical mass of Mdy„«2kpc)=3±lx 10^ c.sc^(/) Mq. 
Adopting a canonical inclination of 30° suggests a dynamical 
mass of yidyn ~10^"Mq. Alternatively, the dynamical mass can 
also be estimated via the velocity dispersion using the relation 
Mdyn=Ca^r /G. The factor C depends on the galaxy mass dis- 
tribution and velocity field, and ranges from C<1 to C >5 (de- 
pending on the mass density profile, velocity anisotropy and rela- 
tive contribution to a from random motion, rotation and assumption 
of spherical or disk-like system; see lErb et alj ( I2006I) for a detailed 
discussion). Adopting C=3.4 (such that a direct compar ison to star- 
forming galaxies at 2=2-3 can be made: lErb et al .l2006h . we derive 
adynamical mass of Md^n ~ 10^" Mq. This dynamical mass is ap- 
proximately 6-lOx smaller than the median Lyman-break galaxy 
mass at z ~3 for which masses have been measured using sim- 



line dynamics (RCS0224-002 arc at z=4.88: ISwinbank et al.l2007h . 
We also note that the inferred dynamical mass is much larger than 
the stellar mass. If the central regions are baryon dominated, this 
may suggest that the gas reservoir makes up >75-90% of the 
baryons. Clearly observations of the cold molecular gas (such as 
the redshifted CO(l-O) or CO(3-2)) would be required to test this. 

Next, we investigate the properties of individual star-forming 
regions within the galaxy. In Fig.|4]we show the extracted the spec- 
tra from the strongest line emitting regions from the galaxy. The 
two largest star-forming regions in Fig[T]have velocity dispersions 
of (Jo=86±15 and a4=80±15 km s~^. Taken with their sizes de- 
rived above we derive masses of star-forming regions of order 6- 
9x10** Mq (we note that in calculating these masses, we have used 
the same value of C as above for simplicity), and star-formation 
rates of SFRo=12±l and SFR4=6±1 M0yr^^ respectively sug- 
gesting that up to half of the total star-formation is occurring within 
two dense regions and outside the nuclear regions. 

Table 1. 

Observed Aperture Photometry for MS 1 358+62arc 



Filter 


mag (AB) 


-B475 


>26.5 


V625 


23.50±0.08 


^775 


22.05±0.05 


2860 


22.12±0.05 


^110 


22.18±0.08 




22.38±0.08 


K 


22.52±0.20 


3.6 /im 


22.90±0.10 


4.5 /im 


23.14±0.I0 


5.8 iim 


23.43±0.25 


8.0 ^tm 


23.50±0.25 



Note. - To convert to intrinsic magnitudes (corrected for the lensing amplifi- 
cation factor 12.5 ± 2. Ox) add 2.7 mags. 



3.4 Comparison with rest-frame UV 

iFranx et al] ( Il997h discuss the morphology and rest-frame UV 
spectral properties of this galaxy in detail. In particular the UV- 
ISM absorption lines show velocity variations on the order of 
200kms~^ along the arc with the SillA1260 line systematically 
blue-shifted with respect to the Lya emission, and a n asymmet- 
ric Ly g emission line with a red tail. As discussed bv lFranx et al.l 
( Il997t) these spectral features are naturally explained by an out- 
flow model, in which the blue side of the Lya line has been 
absorbed by outflowing neutral Hi. The description of the line 
profile is very typical of the emission profiles seen in other 
Lyman-break galaxies. Indeed, velocity offsets between the neb- 
ular emission lines (such as Ha) and UV-ISM emission/absorption 
lines (such as Lya and SillA1260) are now common in high- 
redshift star-forming galaxies, and are usually i nterpreted as ev- 
idence for large-s c ale starburst driven outflow jErb et all l2003l : 
iTeplitz et al] bOOd : IPettini et al.l l2002bl) . Such a model provides 
a good description of the integrated properties of high-redshift 
galaxies and m ore detaile d observations of local starbursts (eg . 



ilar techniques ( eg. lErb et a i] |200d IpHrster Schreiber et al.(l2006l : 
iLaw et al]|2009l) . but comparable to that of the only other star- 
forming galaxy at these redshifts with measured nebular emission 



Pettinietal.ll2002 b: Shaple v et al.ll2003l:lTenorio-Tagle et al.lll999l; 
Heckman et alJl200Q ; iGrimes et alj|2006h . In this model, the Lya 



emission comes from photons emitted from the star-forming re- 
gions. To reach the observer these must pass through part of the 
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foreground (blue shifted) shell. This absorption causes the peak 
emission wavelength to appear redshifted relative to the nebular 
emission lines. In addition, photons may be scattered or emitted 
from the receding shell. Photons that are either created on the inner 
surface of the shell (e.g. by UV irradiation from the starburst), or 
multiply scattered within the receding shell will acquire the mean 
velocity of the outflow and will be seen as red shifted by the ob- 
server (see the discussion of lHansen & oi]|2006h . 

With the systemic velocity of the galaxy measured from the 
[On] emission, we can compare the nebular emission line veloc- 
ity gradient seen in this galaxy with those observed in Lya emis- 
sion and UV-ISM lines, as well as crudely estimating the ener- 
getics of the outflow. Using the systemic redshift measured from 
the [On] emission line doublet {z =4.9296±0.0002) the Lya ap- 
pears systematically redshifted from the systemic redshift by 100- 
200kms~^ across the 4kpc extent of the galaxy (Fig.[5](- More- 
over, the UV-ISM absorption lines are also correspondingly blue 
shifted by ~200±100kms^^. This is typical of high-redshift, star- 
forming galaxies. However, strikingly, the velocity gradient in the 
nebular emission is mirrored in both the Lya and UV-ISM lines. 
The fact that the velocity structure of the outflow follows that seen 
in the nebular emission suggests the outflow has yet to decouple 
from the galaxy. Indeed, assuming that the outflow is no larger than 
the spatial extent of the stellar component of the galaxy (<2kpc), 
then for a constant velocity of 150kms~^ it takes just 15Myr to 
travel 2kpc. This suggests that the young outflows from individual 
Hll regions have not yet merged to form a superwind surrounding 
the galaxy. 

To test whether a young outflow is energetically feasibly, first 
we de-redshift the rest-frame UV spectrum using the systemic 
redshift. In Fig.|4]we show the region around the Lya emission in 
the rest-frame, clearly showing the asymmetric profile. A single 
Gaussian profile fit suggests a velocity offset from the systemic 
redshift of ~200kms~^. However, a much better fit (Ax^ >25) 
is obtained by modeling the emission line with an underlying 
Gaussian profile (fixed in wavelength at I2I5.67A) with an Voigt 
profile absorber in the blue wing. In this fit, the emission line width 
and intensity of the underlying Gaussian profile are allowed to 
vary, as are the velocity, impact parameter (6) and column density 
(n) of the Voigt profile. The best fit model has a Voigt profile 
centered -300±100kms~^ from the systemic with 6 ~ 70kms^^ 
and n ~ 4xl0^^cm~^. Allowing the centroid of the underly- 
ing Gaussian emission line profile to vary, the column density 
increases to n ~ 3xI0^*cm~^ and h ~120kms^^ with an un- 
derlying Gaussian profile with centroid At; = —100 ± 30kms~^ 
from that predicted from the [Oil] emission (we note that the 
improvement in between the two fits is Ax^ < 4 and so 
indistinguishable, although fixing the underlying Gaussian profile 
at the systemic redshift is more physical). We caution that both 
of these column densities are near the flat part of the curve of 
growth causing degeneracies between N^f and b and so come with 
considerable uncertainty. However, we can at least test whether 
such a value is reasonable using the UV absorption lines and 
assuming an ionisation state for the ISM. Using the rest-frame 
UV spectroscopy we measure (rest-frame) equivalent widths 
of Vyo(SiiiAI260.4)=0.30±0.07A VFo(Oil302.I)=I.69±0.07A 
H^o(CiiA1334.5)=0.32±0.07A Wo(SiivAI393.4)=0.24±0.06A 
and V|/o(SilvA1402.8)=0.29±0.06A. Since we can not constrain 
the ionisation state or metallicity of the galaxy (we only have one 
transition from each species), we adopt a simple approach and 
assume that the ionisation and metallicity are comparable to that 
of cB58 aX z = 2.72 which has an abundance of ~ 0.4 solar 



jPettini et alj|2002al) . Assuming log[X/H]=-AM for SillAI260.4 
and an equivalent width of Wo(SillA1260.4)=0.30±0.07A (which 
corresponds to a column density of Nsm ~2x 10^'^ cm~^) we 
derive a hydrogen column density of ~2xI0^^cm~^. Although 
we caution this estimate comes with considerable uncertainty, the 
column density estimated from the Lya line profile, or crudely 
using the UV ISM absorption lines appear, at least, to be consistent 
within a factor ~4x . 

Although there are uncertainties in the column density of the 
outflowing shell, it is interesting to consider the energetics. For a 
normal stellar initial mass f unction, su pernovae provide ~I0*^ erg 
per solar mass of stars (eg.. lBenson e t al. 2003), thus for a galaxy 
with a star-formation rate of SFR=42±7M0 yr"\ ^SxIO^^erg 
are available per year. We estimate the distance of the shell from 
the galaxy by assuming that the linear size of the region covered 
by the swept up shell is no larger than the spatial extent as the 
galaxy (which seems reasonable given the velocity structure of the 
outflow appears to follow the nebular emission line velocity). As 
noted, at a velocity of ^ 150 kms^^ it takes I5Myr to travel out 
to a distance of 2 kpc; within this time, supernovae will provide a 
total energy of Esjvo ~5x 10^^ erg. Clearly there are a number of 
uncertainties in this estimate: not least that the outflow may have 
decelerated from an initial higher velocity which would reduce the 
time-scale. Nevertheless, we can estimate the mass of the outflow 
via: Moutfiow ~ AN hi x tuh /x where Nhi is the column den- 
sity, A is the area of the cone and mn is the atomic mass of Hy- 
drogen. Adopting the value of the observed column density derived 
above (A^'hi ~ x 10^^~^* cm^^) and that the outflowing material 
is uniform over an area of ~10kpc^, the total mass of the out- 
flow is 3 x 10*/a; M0, where x is the neutral Hi fraction. The 
kinetic energy of the outflow is then Ek ~ hmv ~ 

lo^V 

X erg. 

Thus, the outflow is easily energetically feasible even if the neutral 
fraction is as low as 0.01%. The low inferred mass loading and ki- 
netic energy in the shell is in stark contrast to the 2=4.88 galaxy 
behind the lensing cluster RCS0224-002 in which similar observa- 
tions suggest a large s cale, bi-polar outflow surrounds the galaxy 
dSwinbank erZll2007h . In this galaxy, the outflow has a high mass 
loading (comparable or greater than the star-formation rate) and ap- 
pears to be located > 30 kpc from the galaxy, escaping at a speed of 
up to 500 kms^^. The strong contrast between the two galaxies so 
far studied at these early times therefore suggests strong diversity 
in the outflow energetics of young galaxies at high-redshift, clearly 
illustrating the need for more targets and follow-up to test the ubiq- 
uity and impact of outflows at these early times. 

3.5 Comparison of Age and Mass Estimates 

It is useful to compare the various stellar and dynamical mass and 
age estimates, since these reflect the build up of stellar mass in this 
system. First, the stellar mass estimated from the rest-frame UV- 
near-infrared SED of the galaxy is much smaller than the dynam- 
ical mass estimated from the dynamics. Similar results have been 
found in the z ~2-3 population. Erb et al. ( 2006) find ~10% of 
star-forming systems with stellar masses which are approximately 
10% of the dynamical mass. These systems also tend to be the 
youngest; <20Myr (e.g. HDF-BX1397 is an excellent example). 
Thus, although the stellar mass estimate appears small compared 
to the dynamical mass (possibly suggesting a large gas reservoir 
makes up the rest of the baryonic material), the galaxy does not 
appear to be at odds with comparable observations of the galaxy 
population at z ^2-3. 

It is also interesting to note that the age estimated from the 
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Figure 5. Extracted, one-dimensional velocity gradient along the long axis 
of the galaxy (source plane). The velocity gradient of ~180kms~^ ob- 
served in the [Oil] emission is mirrored by the hya and UV ISM lines, 
although they appear systematically offset from the systemic by 150- 
200kms~^, suggesting that the galaxy is surrounded by a young galactic 
scale outflow which has yet to decouple and burst out of the galaxy. Using 
the velocity offsets and spatial extent of the galaxy, we suggest the outflow 
is <Cl5Myr old. 



SED (~15Myr) is comparable to that derived from the instanta- 
neous star-formation rate (at a SFR of 42M0yr~^ it takes just 
~ 1 7 Myr to build a stellar mass ofTxlO^M©). Finally, the dynam- 
ical age of the outflow (~15 Myr) is young, and consistent with a 
scenario in which a starburst wind is in the process of accelerat- 
ing material from the galaxy into the halo and/or IGM. Together, 
these mass and age estimated are all suggestive of a young galaxy 
experiencing its first major epoch of mass assembly. 



3.6 Properties of the Largest Star-Forming Regions 

Owing to the magnification of the galaxy, the source-plane reso- 
lution of the optical and near-infrared imaging we are also able to 
estimate the masses, sizes and luminosities of the brightest star- 
forming (Hll) regions within the galaxy to test whether the phys- 
ical conditions for star-formation at z ^ 5 are very different to 
galaxies today. Using the results from § 3.1 we estimate the bright- 
est star-forming regions within the galaxy have diameters of or- 
der ~300-400 pc from the HST and NIPS observations. Using the 
velocity dispersions measured from the IFU observations this sug- 
gests dynamical masses of ~6-9 x 10* Mq each and mass densities 
of ~3OOOM0 pc~^. In Fig.|6]we show the correlation between size 
and mass of Hll regions (and starburst complexes) in local galax- 
ies. Clearly, the size density of the brightest star-forming regions in 
MS 1358-1-62 arc are consistent with the most massive clusters seen 
in star-forming galaxies. Next we compare the luminosities and 
sizes of the brightest Hll regions to those locally. In Fig.|7]we show 
the corre lation between size and luminosity of local Hll regions 
( Kennic utlll988|:lGonzalez Delgado & Perezll997l : lKennicutt et al.l 
[20031 : ILee et al.(l200 7') and show the position of the two bright- 
est star-forming regions from MS 1358-1-62 arc. As Fig. [7] shows, 
the two brightest Hll regions appear to have star-formation activity 
which is substantially higher (at a fixed size) than those typically 
derived from local star-forming galaxies. 

However, before we interpret this offset in detail, there are two 



caveats which must be considered. We first note that this compar- 
ison relies on a conversion b etween [Oil] ( i n the z — 4.92 arc) 
and Ha locally. However, as Kewlev et al. I ( l2004h point out, the 
conversion between the star-formation rate derived from Ha and 
[On] never varies by more than a factor of two, even for extreme 
ranges in metallicity and reddening. Indeed, targeting the nebular 
emission line r atios of ^^25 lense d galaxies, Richard et al. (2009 in 
prep) (see also Pettini et al. derive a median Ha/[Oll] emis- 

sion line flux ratio of Ha/[Oll]=2.5±0.8 (median redshift of this 
sample z =2.4±0.4), which is slightl y larger, but consist ent with 
the values locally (HQ/rOll]=1.8±0.5: lKewlev et al.ll2004l) . Thus it 
seems unlikely that there is a factor lOOx increase in the [Oll]/Ha 
emission line ratio between z = 2 and z — 5. 

Second, we caution that the isophotes used to extract Hll re- 
gions sizes and luminosities are very different at z =0 compared 
to those at 2 ~ 5 (even accounting for lensing). Hence, to test 
whether we can reliably compare Hll region sizes and luminosities 
at high-redshift to those locally, we make use of the extensive Ha 
narrow-band imaging of local galaxie s, in particular f rom SINGS 
l lKennicuttetai]|2003h and II HUGS jUee et al.ll2007l) . Using 175 
Hq narrow-band images from these samples we test the effect of 
redshifting the z — galaxy images and re-extracting the size- 
luminosity relation. First, we use the the Ha narrow-band imag- 
ing and construct the size-luminosity relation for local galaxies us- 
ing isophotes ranging from 1 x I0~^* to 1 x 10~^^ erg cm~^ s~^ sq- 
arcsec. We find that, as the isophotal flux is increased, the lower- 
luminosity Hll regions fall below the detection threshold, whilst 
only the cores of the higher luminosity Hll regions are identified. 
This results in the local relation remaining in tact, with fewer Hll 
regions identified. In Fig.|7]we show a vector (labeled 'A') which 
shows the typical direction a Hll regions follow as the extraction 
isophote is increased. 

Next, we test the effect of artificially redshifting the imaging 
from a spatial resolution of ~ 10 pc (0.5" at 5 Mpc) to a resolution 
of ~200pc and extract the Hll regions at a fixed isophote. This 
binning has the effect of merging Hll regions, increasing both the 
size and luminosity and in Fig.|7]we show vector 'B' which denotes 
the track which Hll regions follow as the resolution is decreased. 

The simulations show that essentially the size-luminosity re- 
lation does not systematically shift as the images are redshifted, 
and we therefore feel confident that any offset we measure in the 
size-luminosit y relation likely refle cts a real difference their prop- 
erties (see also lMelnick et al.ll2000l .) We note that it is entirely pos- 
sible that the star-forming regions in the z = 5 galaxy represent 
an amalgamation of multiple Hll regions. However, these simple 
simulations suggest that they will remain over-luminous at a fixed 
size. 

Next, we examine the implication of identifying Hll regions 
at z = 4.92 which are significantly more luminous than those 
found locally. At a fixed size of 400 pc, local Hll regions have star- 
formation rates of order SFR=0. IMq yr~^, thus the star-forming 
regions observed in the z —4.92 galaxy appear to be substantially 
more luminous than those in local galaxies. However, at a fixed 
size, offsets of (up to) two orders of magnitude in luminosities are 
observed in the most luminous nearby Hll regions and starburst 
galaxies (eg.|Wil son & Matthews 19951: lEnglish & Freemaij|2003l: 
IVanzi et al.ll2008l: iBastian et al.ll200^ . For example, 30Doradus 
has a star-formation rate of ^O.ISMq yr~^ and a size of ~200pc 
whilst the giant Hll region in the starburst galaxy IIZw40 
has a star-formation rate of ^IMpyr"^ within 400 pc (eg. 
Wil son & Matthews! \l99^ lEnglish & FreemaiJ I2OO3I : IVanzi eFd] 
I2OO8I) . although both of these systems are low metallicity. However, 
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Figure 6. Mass versus size for Hll regions in local galaxies compared to 
the two Hll regions in MS 1358+62 showing that the the two brightest star- 
forming regions are large, massive star clusters, but have comparable den- 
sities to Hll regions found l ocally. The squares denote Milky- Way globu- 
lar clusters (Harr is 1996; Pr yor & Mev lan 1993), inverted open triangles 
represent Ultra-Compac t Dwarfs (Ha^egan et al. 2005; E vstigneeva et alj 
2007tlHilker e t al.l2qq7l) and open triangles are M82 superstar clusters from 
McCrady & Gr ahamI i2007ll. Th e solid triangles denote giant Hll regions 
from iFuentes-Masip et alJi200Cl) . The dashed line shows the predic ted clus- 
ter rad ius as a function of mass (after accounting for mass loss) from lMurrayj 
i2009l) . 

observations of massive star-clusters in NGC 4038/39 (the Anten- 
nae) also indicate young (~5 Myr) clusters with star-formation 
rates much higher t han expected for typical star-forming galaxies 
teastian et al.ll2006h . Thus the Hll regions within this z = 4.92 
galaxy appear to be scaled up versions of the most extreme Hll re- 
gions observed in the local Universe. It is unclear whether this is 
due to the starburst mode, or low metallicity, although it is inter- 
esting to note that the most recent models of galaxy formation sug- 
gest that star-formation at high reds hift may be driv en by strongly 
unstable and fragmenting disks (eg. Dekel et al. 2009) resulting in 
localised starbursts. Thus, the intense star-formation observed in 
this galaxy will significantly effect the surrounding ISM both in 
distributing metals and in adding turbulence. 



4 DISCUSSION & CONCLUSIONS 

After correcting for lensing magnification factor I2.5±2.0x, the 
rest-frame UV/optical morphology appears extended over ~4 kpc, 
but is dominated by approximately five star-forming regions. In 
the source plane the galaxy has an /-band magnitude consis- 
tent with L* LBGs at z~5, and an intrinsic star-formation rate 
of SFR[o//] =42±8M0yr"\ Using rest-frame UV-to-iy-band 
SED we derive stellar mass of M* =(7±2)x lO'^M©, which is 
consist ent with estimates of the stellar masses of LBGs at thi s 
epoch jVerma et al.ll2007l : IStark et alj|2009l ; iMcLure et al1l2009l) . 
The SED is also consistent with being dominated by a young stel- 
lar population of age of 14±7Myr. Indeed, with a star-formation 
rate of ~42M0yr~^, it takes just ~17Myr to build this stellar 
mass. Together, this suggests this galaxy is in the first major epoch 
of mass assembly. 

The source-plane morphology is clumpy with at least five 
star-forming regions, separated by up to 4 kpc in projection. The 



Figure 7. Correlation between size and star-formation rate from Hq for 
Hll regions in local galaxies. The solid points denot e the Hll regions i n 
local star-formin g galaxies from Kennicutt ( 1988); Ke nnicutt et alj j2003h ; 
iLee et alj EoOtI) . We also show the star-formation rate-size relation for 30 
Doradus (the largest, local Hll region in the L MC) as we ll as the giant 
starburst re gion in the nearby galaxy II Zw 40 tVanzi et alj l2008) and the 
Antennae (Bastian et al.ll200g) . The dashed line shows a fit to the Hll re- 
gions in nearby spiral galaxies whilst the dotted hne is the same but shifted 
in luminosity at a fixed size by a factor lOOx. The two Hit regions in 
MS 1358+62 arc are shown by the solid red points where the x error-bars 
denote the range of size estimates from the source-plane IFU and ACS re- 
constructions, and the y-error bar denotes the eiTor in the flux measurement, 
lens modeling u ncertainty and the range of possible [Oii]/Hq emission line 
flux ratios from lKewlev et al] f2004) for extreme ranges in metallicity and 
reddening. The two HlI regions in the z = 4.92 galaxy appear to be scaled 
up versions of the most luminous star-forming regions observed locally, 
even though they are observed when the Universe was less than 1 Gyr old. 



brightest regions are not centrally concentrated, and the galaxy 
shows no signs of a central component. Using the spatially re- 
solved spectroscopy around the [On] emission, we derive a ve- 
locity gradient of Aii=I80±20kms~^ across 4kpc in projection, 
suggesting a dynamical mass of Mdan=3±l X lO'^MQCsc^(i) Mq. 
We also derive a specific star-formation rate (SFR/M*) of 
SSFR=2.4±0.5x 10"** yr"\ consistent with the mos t actively star- 
formi ng, low mass galaxies at these early times ( IFeulner et al.l 
l2005l) . 

Using the source-plane images and [Oil] emission line maps, 
we measure sizes and masses for the two brightest star-forming 
regions, deriving masses of Mc!=6-9x 10* Mq and individual star- 
formation rates of SFR[o//] ^lOM© yr^^. A comparison to lo- 
cal Hll regions/populations shows that these are comparably mas- 
sive and dense as those found in local star-forming galaxies. How- 
ever, the implied star-formation rate individually for these Hll re- 
gions is >I00x brighter that found at a fixed size. Such high star- 
formation rates (at a fixed size) have been observed locally: e.g. 
30 Doradus and II Zw 40 both have star-formation rates which are 
> 100 X larger than expected given their sizes. Thus, the Hll regions 
within this z — 4.92 galaxy appear to be scaled up versions of the 
most extreme regions observed in the local Universe, yet observed 
when the Universe was < 1 Gyr old. Could the increased luminos- 
ity within Hll regions reflect real difference in th e mode of star - 
formation within massive star-forming complexes? iMurravl j2009l) 
suggest that young, massive (>IO^Mo) star-forming clusters are 
optically thick to far-infrared radiation resulting in high gas temper- 
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atures, and hence in higher Jeans masses. Indeed, for a star cluster 
with mass >10* Mq (as observed in this z — 4.92 galaxy) the pre- 
dicted Jeans mass is '^^IIMq. This 'top-heavy' IMF has the effect 
of increasing the fraction of OB stars per Hll regio n and hence in - 
creases the light-to-mass ratio (see the discussion in lMurravll2009h . 
However, since no compelling evidence for evolution in the IMF 
has yet been found at either low (or high) redshift, we view this sce- 
nario as suggestive, at best. Nevertheless, clearly spatially resolved 
spectroscopy of the nebular emission lines on scales comparable 
to Hll regions would provide crucial diagnostics of the physics 
of star-formation in the young Universe. Indeed, this intense star- 
formation will significantly effect the surrounding ISM both in dis- 
tributing metals and in adding turbulence. As such, the mechani- 
cal energy input might explain the higher turbul ent speeds (hence 
large v/g values ) observed in primitive disks (e . g.lStark et alJlOQj ; 
Law et al."2009'; 'Fo rster Schreiber et alj[2009l : lOenzel et alj]200g ; 
Lehnert et al. 2009). 

Finally, we compare the nebular emission line dynamics with 
spatially resolved Lya and rest-frame UV-ISM lines and find that 
this galaxy is surrounded by a galactic scale outflow. The velocity 
gradient observed in the nebular emission is mirrored (but system- 
atically offset) in the Lya emission and UV-ISM absorption across 
~4kpc, suggesting that the outflow is young (<15Myr) and has 
yet to decouple from the star-forming regions and escape from the 
galaxy disk. Although crude, the estimated mass loading in the out- 
flow is very small, and the kinetic energy provided by SNe from the 
observed star-formation rate is easily enough to drive the wind. 

Overall, these observations provide unique insights into the 
distribution of star-formation, dynamics, interaction between star- 
formation and outflow energetics and even properties of the largest 
Hll regions within a young galaxy seen less than 1 Gyr after the 
big-bang on scales of just ~200pc. The combination of the large 
number of high-redshift gravitationally lensed galaxies identified 
by HST and adaptive optics assisted integral field spectrographs on 
eight and ten meter telescopes should finally begin to make such 
critical observational studies common-place, testing the route by 
why early systems assemble their stellar mass, their mode of star- 
formation and how they ultimately develop into galaxies like the 
Milky-Way. 
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